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Abstract. We consider an initial-boundary value problem for the one-dimensional equations of 
compressible isentropic viscous and non-resistive magnetohydrodynamic flows. The global well- 
posedness of strong solutions with general large data is established. Moreover, the vanishing 
resistivity limit is justified and the thickness of magnetic boundary layers is analyzed. The 
proofs of these results are based on a full use of the so-called “effective viscous flux”, the 
material derivative and the structure of the equations. 
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1 Introduction 

Magnetohydrodynamics (MHD) concerns the motion of a conducting fluids in an electromagnetic 
field and has a very wide range of applications in astrophysics, plasma, and so on. Because the 
dynamic motion of the fluids and the magnetic held interacts strongly on each other, both 
the hydrodynamic and electrodynamic effects have to be considered. The three-dimensional 
equations for compressible isentropic magnetohydrodynamic hows, derived from fluid mechanics 
with appropriate modifications to account for electrical forces, read as follows (see & EH): 

p t + div(pu) = 0, 

< (pu) t + div(pu <g) u) + VP = /jAu + (p + A')Vdivu + (V x B) x B, (1.1) 

B* — V x (u x B) = —uV x (V x B), divB = 0 

with x £ A! C M 3 and t > 0. Here, the unknown functions p, u £ R 3 ,P and B £ R 3 are the 
density of fluid, velocity, pressure, and magnetic held, respectively. The viscosity coefficients p 
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and A' satisfy the physical conditions p > 0,3A' + 2p > 0. The constant v > 0 is the resistiv¬ 
ity coefficient which is inversely proportional to the electrical conductivity constant (magnetic 
Reynolds number) and acts as the diffusivity coefficient of the magnetic fields. The pressure 
P(p) is generally determined through the equation of state (the so-called 7 -law): 

P(p) = Ap y with A > 0,7 > 1. (1.2) 

Equations (jl.llh (11.211 describe the macroscopic behavior of electrically conducting com¬ 
pressible (isentropic) fluids in a magnetic field. From Eq. (ll.ljU it is clear that the time rate 
of change of the magnetic held (i.e., Bt) is dominated by both the advection term V x (u x B) 
and the diffusion term z^V x (V x B). However, in many cosmical and geophysical problems 
where the conducting fluid is of extremely high conductivity (ideal conductors), the resistivity 
v is inversely proportional to the electrical conductivity < 7 , and therefore, it is more rational to 
assume that there is no magnetic diffusion (i.e. v = 0) (see, for example, [3, 9]). So, instead of 
(0)3, the induction equation for magnetic held in such cases has the form: 

B t — V x (u x B) = 0, 

which particularly implies that in a highly conducting fluid the magnetic held lines move along 
exactly with the fluid, rather than simply diffusing out. This type of behavior is physically 
expressed as that the magnetic field lines are frozen into the fluid. In effect, the fluid can how 
freely along the magnetic held lines, but any motion of the conducting huid, perpendicular to 
the held lines, carries them with the huid. The “frozen-in” nature of magnetic helds plays very 
important roles and has a very wide range of applications in both astrophysics and nuclear 
fusion theory, where the magnetic Reynolds number R m ~ 1/V is usually very high. A typical 
illustration of the “frozen-in” behavior is the phenomenon of sunspots. For more details of its 
physical background and applications, we refer to mca hi affinal usual. 

Formally, when v = 0, system turns into 

' p t + div(pu) = 0 , 

< (pu)t + div(pu ( 8 > u) + VP = pAu + (p + A')Vdivu + (V x B) x B, (1.3) 

Bi-Vx(uxB) = 0, divB = 0, 

where the pressure P(p) satishes the 7 -law (11,21) . This is often called the compressible isentropic 
viscous and non-resistive MHD equations. 

Because of the complete interaction between the dynamic motion and the magnetic held, 
the strong nonlinear terms and the lack of dissipation mechanism of the magnetic held, many 
physically important and mathematically fundamental problems of system (| 1 ,3(1 are still open. 
For example, to the author’s knowledge, there is no any rigorously mathematical literature on 
the global well-posedness of the initial (boundary) value problem of (11.31) . even that the initial 
data are close to a non-vacuum equilibrium, though the same problem has been successfully 
solved for the compressible Navier-Stokes equations (i.e., B = 0) by Matsumura-Nishida [23]. 

Due to the complicated structure of multi-dimensional equations, instead of (11.11) and (11.31) . 
in this paper we shall consider the simplest one-dimensional equations (cf. [ 12 ]): 

Pt + ( pu) x = 0 , 

< (pu)t + i^pu 2 + P(p) + i b 2 ^j = A u xx , (1.4) 

bt -|- (ub) x — nb xx , 
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and 


Pt + ( pu) x = 0 , 

< {pu) t + (^pu 2 + P{p) + ^b 2 ^j = A u xx , (1.5) 

b t + ( ub) x = 0 , 

where the pressure P(p) obeys the 7 -law ( 11 . 21 ) and A = 2p + A'. 

Indeed, assume that the MHD flow is moving only in the longitudinal direction x and 
uniform in the transverse directions ( y,z ). Then it is easy to derive (11.41) and (11.51) from (11.11) 
and (11.31) . respectively, based on the specific choice of dependent variables: 

p = p(x,t), u = (u(x, t), 0,0), B = (0,0, b(x, t)). 


We mention here that the one-dimensional system for compressible heat-conductive viscous and 
resistive MHD flows in a form similar to that in (11.41) was studied by Kazhikhov-Smagulov H 8 . 
where the global well-posedness of solutions was announced. 

The main purpose of this paper is to show the global well-posedness of strong solutions to 
an initial-boundary value problem of the system (11.51) . to justify the vanishing resistivity limit 
(i.e., v —>• 0) from (II. 4j) to (jl.5j) rigorously, and to study the boundary layer effects as v —>• 0. 
To do so, without loss of generality, we consider an initial-boundary value problem of (11.51) on 
a bounded spatial-domain D = ( 0 , 1 ) with the following initial and boundary data: 


(p,u,b)(x, 0) = (p 0 ,u 0 ,b 0 )(x), x G [0,1], 
u(0,t) = u(l,t) = 0, t> 0. 


( 1 . 6 ) 


Our first result concerns the global well-posedness of strong solutions to the initial-boundary 
value problem (11.51) . (11.61) . 

Theorem 1.1 Assume that the initial data ( po,uo,bo ) given in \1.6 \ ) satisfies 

inf p 0 (x) > 0, (po,b 0 ) G H 1 , u 0 G Hq D H 2 . (1.7) 


Then for any 0 < T < 00 , there exists a global unique strong solution (p, u, b) to the initial¬ 
boundary value problem U.5\) . \1.6\ ) on [0,1] x [0, T), satisfying 

0 < C~ x < p(x, t) < C < 00 , V (x, t ) € [0,1] x [0, T) 

for some positive constant C, and 

(p, b) € T°°( 0, T; H 1 ), (p t , b t ) € L°°( 0, T; L 2 ), 
u G L°°(0,T; Hq D H 2 ), u t G L°°(0, T; L 2 ) D T 2 (0, T; H 1 ). 

Theorem 1 1.1 1 is the first result on the global well-posedness theory of the non-resistive MHD 
equations with large data. It particularly implies that the solutions of the viscous and non- 
resistive MHD problem (11.51) . (11.61) will not develop vacuum and mass concentration in finite 
time provided the initial data are bounded and smooth, and do not contain vacuum. 

It is worth pointing out that the equations m look similar to a compressible model 
for gas and liquid two-phase fluids (see, for example, PSIE]). However, to prove the global 
well-posedness of the two-phase model, it is technically assumed in [291 [ 6 ] that the proportion 


( 1 . 8 ) 

(1.9) 
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between the mass of gas and liquid must be bounded, in analogy to the assumptions p,b > 0 
and 0 < b/p < oo for ill.51) . Of course, this is somewhat reasonable for the two-phase model, 
but not physical and realistic in magnetohydrodynamics. 

The second purpose of this paper is to justify the vanishing resistivity limit from (11.41) to 
(11.51) rigorously, as the resistivity coefficient v —>• 0. More precisely, 

Theorem 1.2 (i) Suppose that system o is equipped with initial and boundary data: 


f {p,u,b)(x, 0 ) = (p 0 ,u 0 ,b 0 ){x), 

x € [ 0 , 1 ], 

( 1 . 10 ) 

|u( 0 , t) = «(1, t) = 0 , 6 ( 0 , t) = 

blit), 6 ( 1 , 6 ) = b 2 (t), t> 0 . 

For any 0 < T < oo, assume that 



inf Po{x) > 0, (p 0 , b 0 ) € H l , 

v<x<l 

( 6 1 : 6 2 ) € C 1 ([0, T)), u 0 G Hq. 

( 1 . 11 ) 


Then the initial-boundary value problem hi. 1(A) has a global unique strong solution (p, u, b) 

on [0,1] x [0, T), satisfying 

0 < C^ 1 < p(x, t) < C < oo, V (x, t ) € [0,1] x [0, T) (1.12) 

and 

sup C||«x ||£2 + + V 1/2 \\px\\l 2 + v 1/2 \\b x \\ 2 L 2 ) (t) 

0 <t<T v ' 

+ J (ll/° 1 ,/ 2 ^lll 2 + v z l 2 \\bxx\\ 2 L2 + dt < C, (1-13) 

where C is a positive constant independent of v. 

(ii) Assume that {p v ,u v ,b v ) and ( p,u,b), defined on [0,1] x [0, T), are the solutions of the 
problems | l.f\ ), hl,10 \ ) and 11.6]) . respectively. Then, 

[ {p v , u v , b v ) —> (p, u, b) strongly in L°°(0,T;L 2 ), 
uf —» u x strongly in L 2 (0,T;L 2 ), 

and moreover, there exists a positive constant C, independent of u, such that for v S (0,1), 

sup {\\p v - pf L 2 + W - h\\ 2 L2 + \\u v - u\\ 2 L2 ) it) + [ \\{u v - u) x \\ 2 L 2 dt <Cu 1/2 . (1.14) 

0 <t<T JO 

The third and main result of this paper is to study the effects of magnetic boundary layer 
as the resistivity coefficient v —>• 0. In fact, it is obvious that when the resistivity coefficient 
goes to zero, the parabolic equation (11.41) -^ turns into the hyperbolic equation (J1.5IU . and the 
boundaries become characteristic due to the non-slip boundary conditions ti| x =o,i = 0- Thus, 
by the classical theory in [ 21 ], one has to drop the boundary conditions of the magnetic field 
in ( 11 . 101 ) (cp. (11.61) 1. and consequently, because of the disparity of boundary conditions, we 
cannot expect that as v —>• 0, the solution of the problem (11.41) . (11.101) will tend to the one of 
the problem (11.51) . (11.61) uniformly up to the boundaries x = 0,1. In other words, a (magnetic) 
boundary layer appears near the boundary. 

Similar to the relations among the Euler, Navier-Stokes and Prandtl equations (see, for 
example, | 20 , .25] GSj), it is expected that as v —» 0 , the solution of the problem (11.4p . ( 11 . 101 ) 
converges uniformly to the solution of the problem (11.51) . (11.61) away from the boundaries, while 
there is a sharp change of gradient near the boundary. Inspired by this, we introduce the concept 
of magnetic boundary-layer thickness (MBL-thickness) as follows. 


4 














Definition 1.1 A non-negative function 5{u) is called a MBL-thickness of the problem 
hi. Ill) with vanishing resistivity limit, if 6(v) | 0 as v f 0, and 

lim II ( P u ~P,u u - u, b u - 6 ) || i oc ( 0 i r ; c(n, (l/) )) = °> ( L15 ) 

\\(p u - p,u u - u,b u - b)\\ LOO{0J ,. c{U)) >0 (1.16) 

with = {x € D | 6 < x < 1 — 5}. Here, {p u ,u v ,b u ) and ( p,u,b) are the solutions of the 
problems hi. 1C\ ) and hl.5\ ) . hl.6\) . respectively. 

The concept of boundary-layer thickness (BL-thickness) has been introduced in | T 0 l fTTTI 
in a similar manner for the one-dimensional cylindrical compressible Navier-Stokes equations 
with vanishing shear viscosity limit. The BL-thickness for the scalar conservation laws and 
the 2D Boussinesq equations with vanishing diffusivity limit was also studied in m and [ 15 ] . 
respectively. It is worth mentioning that Definition 11.11 does not determine the MBL-thickness 
uniquely, since any function satisfying 6 i {u) > 5(is) and 5\ (u) | 0 as u j. 0, is also a MBL- 

thickness. Thus, there should exist a minimal MBL-thickness 5*(y) which may be considered as 
the true MBL-thickness. 

In this paper, we shall prove that a function 5 n {y) = iA/ 2 ^ 1 / n with n > 2 is a MBL-thickness 
in the sense of Definition 11.11 This is somewhat in agreement with the famous Stokes-Blasius 
law in the laminar boundary layer theory (see m), since liminfn^oo S n (u) = v 1 ! 2 . In order to 
simplify the analysis of MBL-thickness, similarly to mm, we shall focus on the special case 
of vanishing initial data. More precisely, we shall prove 

Theorem 1.3 (i) Assume that the viscous and non-resistive system \1.5\ ) is equipped with initial 
and boundary data: 

(p, u, b)\ t= Q = (p, 0,0) with ~p = Const. > 0, and u\ x= op = 0. (1-17) 

Then the problem hi. ,51) . jl.l r /\ ) has only a trivial solution ( p,u,b ) = (p, 0,0). 

(ii) Assume that ( p,u,b ) is the solution of the viscous and resistive system (Z3P with the 
following initial and boundary data: 


(p,u,b)\ t=0 = (p, 0 , 0 ), u\ x=0}1 = 0 , 6 ( 0 ,f) = 6 i(f), 6 ( 1 , t) = b 2 {t), (1.18) 

where ~p = Const. > 0 and b\(t), 62 (f) are the same as in |1.17| ) and hi. 101) . respectively. Then 
any function 5(v) > 0, satisfying 


5{v) —>• 0 and . —>• 00 , as v —>• 0 , 
v v t 1 

is a MBL-thickness in the sense of Definition ) 1. 1\ such that 


0 


(IIP — /^ll 


L-(0,T;C(f2 5M )) 


+ 


lL-(0,T;C(n 5M )) 


= 0 


and 


“ (\\P Plli~( 0 ,T;C(n)) + 


L°°(0,T;C(f2))y > 


provided the boundary data 61 (i), 62 (t) are not identically zero. 


(1.19) 


( 1 . 20 ) 

( 1 . 21 ) 
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In view of Definition 11.11 we notice from (11.2011 and (11.211) that there is no boundary layer 
effect on the velocity. This is mainly due to the smooth mechanism of the viscosity term A u xx 
with A > 0, and it indeed holds that (see (|4.4D . (14. 5 p below) 

H' U llL oo (0,T;C(n)) “ aS V 

The proofs of Theorems II. 1111.31 will be given respectively in Sections 1, 2 and 3, based 
on the global (uniform-in-^) a-priori estimates of the solutions. However, the lack of smooth 
mechanism of the magnetic field, the strong nonlinearities, and the interaction between dynamic 
motion and magnetic field will cause some serious difficulties. To circumvent these difficulties, 
motivated by the study of multi-dimensional Navier-Stokes/MHD equations (see, for example, 
[21 £21 EH), we shall make a full use of the so-called “effective viscous flux” F and the material 
derivative u: 

b 2 

F = A u x — P(p) —— and u = ut + uu x . (1.22) 

It turns out that the “effective viscous flux” F possesses more regularities than the velocity. 

We begin the proofs with the elementary energy estimates and the upper boundedness of the 
density (see Lemmas 12.1112.2113.21 and 13.31) . Next, by using the special mathematical structure 
of F. especially, using the non-negativity of b 2 and the fact that F x = pu due to (I1.4IL (or 
ccsy, we can improve the integrability of the magnetic field and obtain the desired bounds 
of \\u x \\ L °°(o, T[L 2 ) and \\p 1,,2 u\\l 2 (o,T;L 2 ) ( see Lemmas I2.3H2.5I and !3.4p . With these estimates at 
hand, we then can show the boundedness of the magnetic field and the lower boundedness of 
the density as well (see Lemmas 12.61 and 13.51) . Thus, noting that F x = pu , by (jl .221) we deduce 

Il'Mxllz, 00 < c (||F|| L oo \\P(p)\\ L oo + 11b"11Z/°°) < C (1 + ||F|| l ~) 

<C( 1 + ||F|| L 2 + \\F x \\ L 2 ) < C (l + \\p 1/2 u\\^) € L 2 ( 0 ,T), (1.23) 

which plays a very important role in the entire analysis of this paper. Indeed, as an immediate 
result of (11.231) and the blowup criterion in [28] . one can easily obtain the global existence of 
strong solutions to (11.51) . (11.61) (see Section 2). 

The justification of vanishing resistivity limit and the study of magnetic boundary layers are 
more difficult and need some more delicate estimates. Indeed, due to the presence of boundary 
layer effects, the global uniform-in-i/ estimates for the solution (p,u,b) of the problem (II.4p . 
(I1.10D (also (11.41) . (|1.18j) ) are much fewer than those for the solutions of the problem (11.51) . (11.61) . 
For example, it is very difficult to obtain the global uniform (in u) L 2 -bounds of the derivatives 
of the density and magnetic field for the problems (11.41) . (11.101) and (11.41) . (11.181) . compared with 
the uniform bounds stated in Lemma [2771 for the problem (11.51) . (11.61) . Instead of (12.271) . noting 
that 

t^xx — A (Fx T P(p)x T bbzc') , 

using (11.231) . the estimates obtained and subtle boundary analysis, we have (cf. Lemma [3761) 

u l/2 sup (\\px(t)\\ 2 L 2 + \\b x (t)\\ 2 L 2) +V 3/2 [ \\b xx \\ 2 L2 dt < C. (1.24) 

0 <t<T Jo 

This suffices to prove the vanishing resistivity limit and to obtain the convergence rates given 
in Theorem [L21 (see Section 3). 

As for the analysis of boundary layer effects, the weighted L 1 -method used in PIS 25] 
seems difficult to apply here, due to the strong interactions of dynamic motion and magnetic 


6 


























field. Instead of the L 1 -method, we shall make use of the weighted L 2 -method to analyze the 
thickness of magnetic boundary layers. To do this, we first utilize the special initial data in 

(11.181) . (11.141) and ([1.241) to improve the convergence estimate of the velocity (cf. Lemma 14.11) : 

sup \\u x (t )\\ 2 L 2 + [ ||/ 0 1 /, 2 u||^ 2 dt < CV 1 / 2 . (1.25) 

o <t<T Jo 

This particularly implies that ||w(t)|| c ^ —>• 0 as v —» 0 for t € [0, T), and hence, there is no 
boundary layer effects on the velocity. Using (11.221) . (II.240 . ([1 ,25[) and the special initial data in 

(11.181) again, we succeed in deriving the weighted (interior) T 2 -estimate (cf. Lemma 14.21) : 

sup f £(x) (\p x \ 2 + \b x \ 2 ) (x, t)dx < Cv 1 ^ 2 with £(x) = x 2 (l — x) 2 , (1-26) 

o <t<T Jo 

which, together with (11.141) and Sobolev’s inequality, proves Theorem 11.31 (see Section 4). 


2 Global well-posedness of (11.51) . (11.61) 


The local existence of regular solutions with smooth initial data can be shown by the standard 
method based on the Banach theorem and the contractivity of the operator defined by the 
linearized the problem on a small time interval. The global existence of solutions will be obtained 
by the method of extending local solutions with respect to time based on the global a-priori 
estimates. 

For this purpose, let ( p,u,b ) be a smooth solution of ( 11 . 51 ) . ( 11 . 61 ) . We will establish the 
necessary global a-priori estimates of ( p , u, b ) defined on [0,1] x [0, T ) for any fixed T > 0. For 
simplicity, we denote by C and Ci (i = 1 , 2 ,...) generic positive constants which may depend 
on A, A, 7 , the initial norms of (po,uo,bo) and T, and may change from line to line. 

We begin with the conservations of mass and momentum. 


Lemma 2.1 Let (ppu,b) be a smooth solution of 11.51) . 11.61) on [0,1] x [0,T). Then for any 
0 <t<T, 

p(x,t)dx = / po(x)dx < oo, (2-1) 

Jo 

and 



j ^pu 2 + ^b 2 +—^-jpA (x,t)dx + \ J \\u x \\ 2 L2 ds 

= So (\ poU ° + ^ + t^T^) ^ dx ' 


( 2 . 2 ) 


The upper boundedness of the density can be deduced in a similar manner to that in [7]. 
Lemma 2.2 Let ( p,u,b ) be a smooth solution of 11.51) . 11.61) on [0,1] x [0,T). Then, 


0 < p(x,t) < C, V (x,t) € [0,1] x [ 0 ,T). 


(2.3) 
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Proof. The non-negativity of the density (i.e., p > 0) readily follows from the method of 
characteristics and the fact that po > 0. In the next, for completeness we sketch the proof of 
the upper bound below. Define 


pt / ^2 \ /*x 

ip(x,t) = J \\u x -pu 2 - P{p) -—J (x,s)ds +J (po«o)(£) d £- 
Clearly, it follows from () 1.5 [) i and (11.51) 9 that 

b 2 f x 

i>x = pu, ift = A u x - pu 2 - P(p) - —, if\ t =o = / (pouo)(£)d£, 

z Jo 

and hence, by Lemma 12.11 one has 

/ if(x,t)da 

Jo 


(2.4) 


(2.5) 


< C and \\iI>x\\l°°(p,t-,l1) < C, 


which particularly yields 


L 00 (0,T;L 00 ) < 


/ if(x, t)da 

Jo 


+ II Vh: IIl oo (0,T;L 1 ) < C. 


( 2 . 6 ) 


Let D t = dt + ud x denote the material derivative and set 

if(x,t) 


(2.7) 


4>(x,t) = exp 
Then, by straightforward calculations we have 

D t (p&) = d t {p<&) + ud x {p§) = ~ (p{p) + ^J p® < 0, 

and consequently, 

l|(^)(i)IU»<||(^)(0)|U=o<C', 

which, combined with the fact that C~ l < 4?(x,t) < C for all (x,t) G [0,1] x [0, T) due to (12.61) . 
immediately leads to (12.31) . □ 

Due to lack of dissipation mechanism of the magnetic field and the strong coupling of dy¬ 
namic motion and magnetic field, the lower boundedness of the density, which relies strongly on 
the L°°-norm of magnetic field b , is more difficult to achieve, compared with the upper bound. To 
circumvent the difficulties, motivated by the mathematical theory of multi-dimensional Navier- 
Stokes/MHD equations (see, for example, [ 8 l [22[ [21]), we introduce the so-called “effective 
viscous flux” F, which possesses more regularities than the velocity. Define 

F(x,t) = ^\u x — P(p) — (x,t) and u(x,t) = (u t + uu x )(x,t), ( 2 . 8 ) 

where F and “ ’ ” are the “effective viscous flux” and the material derivative, respectively. 
Thanks to (j 1.5 j) i and (11.51) 9. it is easy to see that 

F x = p(u t + uu x ) = pu. (2.9) 


The quantities of “effective viscous flux” F and material derivative u will play an important role 
in the entire analysis, particularly, in controlling the first-order derivatives of the solutions and 
studying the vanishing resistivity limit and the boundary layer effects. 

In order to estimate the first-order derivative of the velocity, we first need the following 
preliminary lemma. 







Lemma 2.3 Let F and ” be the same as in \2.8\) . Then for any 0 < t < T, 

x d '\h + ilp 1/2 »fe < d h 1 1 C -■ ■>- • ^ '""“ 3 


2 dt 


dt 


P{p) + ~2 ) u xdx + C\ (Hollis + lifeline) + C. 


Proof. Multiplying ( 11 . 51) 9 by u in L 2 and integrating by parts, we deduce 

X d C D C 1 

/ u 2 dx + / pu 2 dx = — / P(p) x (u t + uu x )dx 
1 Jo Jo Jo 

1 C A l 1 

-- / {b 2 ) x (u t + uu x )dx - - / v? x dx. 

1 Jo 1 Jo 

Due to (USD i, it holds that 

P(p)t + uP(p) x + 7 P(p)u x = 0 with P(p) = Ap 1 , 
and hence, the first term on the right-hand side of (12.111) can be estimated as follows. 


[ P{p)x(u t 

Jo 

r p(p) 

Jo 


+ uu x )dx = / ( P(p)u xt - P(p) x uu x ) dx 


/o 


d_ 
dt 

d rl 


»i 


u x dx - I ( P{p) t u x + P{p) x uu x ) dx 
l 


10 


dt 


r o 


P(p)u x dx + 7 / P(p)u 2 dx. 

Jo 

Similarly, due to (I1.5ID . one has 

(& 2 )t + u(b 2 ) x + 2Fu x = 0 , 


so that, 


1 r 1 If 1 

- / (b 2 ) x (u t + uu x )dx = - (b 2 u xt - ( b 2 ) x uu x ) dx 

2 Jo 1 Jo 

f b 2 u x dx-- f ((b 2 ) t u x + (b 2 ) x uu x ) dx 

Jo * Jo 


1 d 

2dt 

1 d_ 

2dt 


f b 2 u x dx + f b z ufdx. 

Jo Jo 


2 „. 2 , 


By virtue of (12.31) . we infer from (12.81) that 

\u x \ 3 <C(1 + \F\ 3 + b 6 ), V (x,t) € [0,1] x [0,T). 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 


(2.15) 


(2.16) 


Therefore, substituting (|2.13l) . (12.151) into (12.IIP , and using (12.31) . (I2.16|) and Cauchy-Schwarz’s 
inequality, we immediately obtain ( 12 . 101 ) . □ 

The next lemma is concerned with the higher integrability of the magnetic field. 

Lemma 2.4 Let F and “ ’ ” be the same ones defined in \2.8\) . Then for any 0 < t < T, 

d 


dt 


I l 4 + II^IIl 6 — c + c\\f\\ l3 . 


(2.17) 
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Proof. Indeed, multiplying (11.5|h by 46 3 and integrating by parts, we get from (12.81) that 

■f 


0 = — f b 4 dx + 3 / u x b q dx 
dt Jr\ 


= 4 - [ b 4 dx + v 


dt 


A 


b 2 


F + P(p) + j ) b 4 dx 


+ + < F + p W)'' 4 ^ 


which, combined with (12.31) and Cauchy-Schwarz’s inequality, leads to (J2.17I) . 

Combining (12.101) with (12.171) . we easily obtain 
Lemma 2.5 Let ( p,u,b ) 6 e a smooth solution of U.5\) . \1.6\) on [0,1] x [0, T). Then, 

(ll+rWIll 2 + ||6(f)||| 4 ) + [ 

o <t<T 

and moreover, 


□ 


sup (|M 6)||| 2 + || 6 (t)||| 4 ) + / (\\p l,2 u\\h + \\b\\%)dt<C, (2.18) 

)<t<T Jo v ' 


0 <t<T 

where F and ” are defined in (ED 


sup \\F(t )\\ 2 L 2 + / \\F x \\ 2 L2 dt < C, 
<t<T Jo 


(2.19) 


Proof. Multiplying (12.171) by a (large) number /<" > max{(2A) 1 , 2C{\ and adding the resulting 
inequality to (I2.10p . we deduce 


d 

dt 


J (^u 2 c + {K—( 4A) 4 ) b 4 ^ dx + J [pu 2 + (K — C\)b 6 ) dx 
< | f Firing - \± f (yx«, - +c l' |F|^ + C, 


which, integrated over ( 0 , t), gives 

( 11 ^( 6 ) 11^2 + I| 6 (t)||| 4 ) + (||p 1 / 2 'it ||^2 + llftllie) ds < C + C 


[Fills ds, 


( 2 . 20 ) 


since integration of the second term on the right-hand side is non-negative and it follows from 
(12.31) and Cauchy-Schwarz’s inequality that 

A,, ,,Q ,-vllT-./ x I I 2 A, 


— g ll'+z ||| 2 + C||P(p )||| 2 < — ||+r||| 2 + C. 


P(p)u x dx 
3 

Thanks to Lemmas o and 12.21 we infer from (12.81) and (|2.9I) that 

|F || L 1 < C (HuxIIl 1 + \\p\\li + II 6 III 2 ) < C (1 + H^xIIl 2 ) , 

|F||| 2 < C (l + \\u x \\ 2 l2 + ||6||| 4 ), ||Fx||| 2 < C , 2||p 1/2 'u||| 2 - 

This, together with Sobolev’s and Cauchy-Schwarz’s inequalities, yields (0 < e < 1) 
11 F 111 3 < ||F|| L i ||F || 2 cxj < <7(1 + ||ttx||L 2 ) (||F || 2 2 + 11 F 11 jji 11 F x 11 ^2 ) 


( 2 . 21 ) 
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< C*2 1 e||-F 3; ||£2 + Ce 1 (1 + 11^11^2)11-^11^2 

— 4p lf ^\\h Cz 1 (1 + 11 11 ^2 ) (1 + ||Wj;||^2 + ll&ll^u) • (2.22) 

Therefore, substituting (12.221) into (12.201) . and choosing e > 0 sufficiently small, by Gron- 
wall’s inequality we obtain (12.181) . since (12.21) implies that ||itx(i)|||2 € L 1 (0,T). As a result, 
(12.191) follows from (12.211) and ([2.181) . □ 

Next, by making a full use of the mathematical structure of “effective viscous flux” F 
again, we can prove the global boundedness of the magnetic field, and consequently, the lower 
boundedness of the density can be derived in a manner similar to that used in the derivation of 
the upper boundedness in Lemma 12.21 

Lemma 2.6 Let (p,u,b) be a smooth solution of hl.5\) . hi. (A) on [0,1] x [0,T). Then, 

sup (||&(£)||l°o + ||p —1 (f)|| l°°) < C. (2.23) 

0 <t<T 

Proof. Multiplying (I1.5IL by 2n6 2n_1 with 1 < n E N, integrating by parts over (0,1), and 
recalling the non-negativity of P(p) and b 2 , we obtain 

— f b 2n (x, t)dx = — (2n — 1) f u x b 2n dx 
dt Jo Jo 

=- 5 x i r( F+p(rt+ 9 t2 ” <fa 

2r7 — 1 f 1 

<- — Fb 2n dx<Cn\\F\\ Lao \\b\\ 2 ff n . (2.24) 

A Jo 

It follows from (12.31) . (12.91) . (12.181) . (12.191) and Sobolev’s inequality that 

\\F\\l°° < C {\\F \\ 2 L 2 + 11^11^11^11^) <<7(1 + ||p 1 / 2 u|| L 2 ) g L 2 (0,T), (2.25) 

and hence, by direct calculations we deduce from ()2.24|) and (12.251) that for V t € [0, T], 

115(t)11 2 jffn < C 2n exp / Cn J ||.F||£°°dt| < C 2n exp {Cn} , 

where C is a positive constant independent of n. Thus, if we raise to the power l/(2n) to both 
sides and let n —>• 00, then we get 

116(4)11^00 < C, V t £ [0, T). (2.26) 

The lower boundedness of the density can be shown in the same way as in Lemma 12.21 Let 
ip and be the same functions defined in (12.41) and (12.71) . respectively. Analogously to the proof 
of Lemma 12.21 we find 

Dt (yi) =+ V) ^ with 

so that, using (12.31) . (12.61) and (12.261) . we have 

||(p4 , r 1 (t)|| Z/00 < ||(p4>) _1 (0)|| ioo exp|i^ (\\P(p)\\l°° + ^||b|||oc^ dij < C, 

which, combined with (12.6|) again, yields a desired lower bound of the density. □ 

The next lemma concerns the first-order derivatives of the density and magnetic field. 
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Lemma 2.7 Let ( p,u,b ) be a smooth solution of \1.5\) . 47.61) on [0,1] x [0, T). Then, 


SU P (IIAeIIl 2 + HMl 2 “1“ llptII l 2 + INIl 2 + |M|z,°°) ( t ) < C. (2.27) 

o <t<T 

Proof. First, in view of (12.31) . (12.81) . (|2.19l) . (12.231) and Sobolev’s inequality, we have 

ll^hcllz, 00 < c (||f|| L oo + \\p(p)\\ LOO + 

< c (1 + ||F|| L oo) < c (1 + \\F x \\ L 2 ) (2.28) 


ll^xxllz^ < C (\\F x \\ L 2 + \\P{p)x\\l 2 + II(^ 2 )x||l 2 ) 

< c (II^eIIl 2 + IIpxIIl 2 + II^IIl 2 ) ■ (2.29) 

Next, if we differentiate (U.5D i . (|1.5|U with respect to x, multiply the resulting equations by 
p x ,b x in L 2 respectively, and integrate by parts, then we deduce from (12.31) . (|2.23|) . (12.281) and 
(12.29[) that 


fif {\\px\\ 2 L 2 + IIMi 2 ) 

< C'II^Hl 00 (||px|||2 +11 11 2 2) + c|| ^xx ||l 2 (IIpxIIl 2 + HM L 2 ) 

< C (1 + H^xlU 2 ) (11 Px 11^2 + H&xlli 2 ) 

+c{\\f x \\ L 2 + \\p x \\ L 2 + H&xIIl 2 ) (IIpxIIl 2 + II^xIIl 2 ) 

< C (1 + ||F x || l2 ) (l + \\p x \\ 2 L 2 + 11 b x ||^ 2 ) , (2.30) 

which, combined with (12.191) and Gronwall’s inequality, shows 

sup {WPx.mh + \\b x (t)\\l 2 ) < C. 


Moreover, it readily follows from (1 1.5 D i and (11.51) 2 that ||/t?£(*)Hz, 2 , HMOIIl 2 are bounded on [0, T), 
since one easily deduces from (12.21) . (12.31) . (12.181) . (12.231) and Sobolev’s inequality that 

IK*)IIl« < C (\\u{t)\\ 2 L2 + \\u(t)\\ L 2 \\u x (t)\\ L 2 ) <C, V t G [0,T). 

The proof of Lemma 12.71 is thus finished. □ 

Based on Lemmas l2.5H2.7l we can easily derive higher-order estimates for the velocity. 

Lemma 2.8 Let ( p,u,b ) be a smooth solution of 111.5]) . 47.61) on [0,1] x [0, T). Then, 

(ll^xxlli 2 + H^xIIloo + ll^tlli 2 ) dt < C , (2-31) 

and moreover, 

sup (\\u xx (t )\\ 2 L 2 + IKWIli 2 ) + / \\u xt \\ 2 L2 dt<C. (2.32) 

0 <t<T JO 
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Proof. First, it follows from (12.31) . (12.81) . (12.91) . (12.181) . (12.191) . (12.231) and (|2.27[> that 


[ (ll^xxllia + ||?+|l| 2 ) dt < C f (\\Fx\\ 2 L2 + \\P(p)x\\ 2 l2 + \\{b 2 ) x \\ 2 L2 ) dt 

Jo Jo 


u | \j^2 d - ||n| 


L°° II u x || 1 , 2 ) 


< 


ref 

c + c f o {\\F x \\ 2 L2 + \\ P 1/2 u\\ 2 L ,)dt<c, 


which, together with ()2.18|) and Sobolev’s inequality, yields that llri^H^oo € L 1 (0,T). 
To prove (12.321) . we differentiate (11.51) 9 with respect to t to get 

/ 6 2 \ 

putt + puu xt - A u xxt + \ P(p) + — j = -pt(ut + uu x ) - pu t u x , 

which, multiplied by ut in L 2 and integrated by parts, results in 

1 d rl 


\u^s: 


pufdx + A I ufdx 
o 

1 r / 5 2 \ 

P(p) + yj u xt - pt(v% + uu x u t ) - pufu x 


/o L 


dx 


< C 


A(N + N) 

Jo 


U xt | + \u t \{\pt\\u t \ + \pt\\u x \ + \u x \\ut\)} dx 


< c (||pt||L 2 + IIMl 2 ) ll+rtllz, 2 + C (||lh||| 2 + 11 Pt. 11 < \,2 + ||«x|| 2 2 ) 

< C\\u xt \\ L 2 + C (\\u t \\ L 2 + 11 U x t 11X /2 ) (l + ||«i|| 2 2 ) 

< -j\\ U xt\\ 2 L 2 +(7(1+ ||ttt||| 2 ) , 

where the previous lemmas and Cauchy-Schwarz’s inequality have been used. Thus, it follows 
from (I2.23p . (12.311) and Gronwall’s inequality that 


f T 

sup \\u t (t)\\ 2 L 2 + / \\u x t\\ 2 L2 dt < C. 
0 <t<T Jo 


As a consequence, using (12.31) . (12.181) . (12.231) and (12.271) . we infer from (11.51) 9 that ||'Wxx(^)||l 2 is 
bounded on [0,T). The proof of Lemma 12.81 is thus finished. □ 

Based on the local existence result and the global a-priori estimates established in Lemmas 
I2.1H2.81 we can prove Theorem II. 11 

Proof of Theorem 11.11 As aforementioned, the local-in-time solutions can be obtained via 
the standard fixed point theorem. Thus, based on the global a-priori estimates established in 
Lemmas 12.1112.81 we can extend the local solutions globally in time on [0, T ) for any 0 < T < oo. 
This proves the global existence of strong solutions. The uniqueness of strong solutions can be 
easily shown by using the standard L 2 -method, and the details are omitted here for simplicity. 
The proof of Theorem 11.11 is therefore complete. □ 
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3 The vanishing resistivity limit 


The global existence of strong solutions to (|1.4jl , (11.101) stated in the first part of Theorem 11.21 
can be established in a manner similar to (indeed much easier) that used in [5) [30] by combining 
the standard local existence result and the global a-priori estimates. Thus, in this section we 
only focus on the derivation of the uniform bounds stated in (11.121) and (11.131) , which will suffice 
in the study of the vanishing resistivity limit in the second part of Theorem ll.21 To do this, we 
assume that (p, u , b ) is a smooth solution of (II. 4|) . ([l.lOp defined on [0,1] x [0, T). For simplicity, 
we also denote by C a generic positive constant, which may depend on A, A,^y,T, the norms of 
the initial data (po,uo,bo)(x) and boundary data (bi,b 2 )(t), but is independent of v. 

First, to deal with the boundary effects, we need the following formulas of b x on the bound¬ 
aries x = 0,1. 


Lemma 3.1 Let ( p,u,b ) be a smooth solution of jjl.ff) , 111. 1(A) on [0,1] x [0, T). Then, 


vb x (0,t) = v(b 2 (t) - bi(t)) - d t 
vb x (l,t) = v(b 2 (t) - bi{t)) +d t 


1 / px 


fe(£, ) dx 


i o vo 

r 1 / 


fe(£, t)d£ ) dx 


' 0 \J x 


/ ( ub)(x,t)dx , 
Jo 


— / ( ub)(x,t)dx. 


(3.1) 

(3.2) 


Proof. First, thanks to the non-slip boundary condition rt| x= o = 0, we deduce after integrating 
Da over (0, x) that 

vb x (0,t) = ub x - d t [ b(£,t)d£-ub, (3.3) 

Jo 

and hence, by integrating ([3.3D with respect to x over (0,1), we obtain ()3.1D . 

Similarly, integrating (jl.4|U over (x, 1) and using the non-slip boundary condition u\ x= \ = 0 
again, we find 

vb x (l,t) = nb x + d t [ b(£,t)d£-ub. (3.4) 

Jo 

Thus, an integration of (13.4|) with respect to x over (0,1) immediately leads to (13.21) . □ 

In terms of the boundary formulas stated in Lemma [3.11 one can easily derive the elementary 
energy estimates similar to the ones in Lemma 12.11 


Lemma 3.2 Let ( p,u,b ) be a smooth solution of il,10\) on [0,1] x [0, T). Then for any 

0 <t<T, 

/ p(x,t)dx = / po(x)dx, (3-5) 

Jo Jo 

and 

J {)^P u2 + \b 2 + (x,t)dx T J (\\\u x \\ 2 L 2+v\\b x \\ 2 L2 )ds <C. (3.6) 


Proof. As a consequence of conservation of mass, one obtains (13.51) . To prove (I3.6jl . multiplying 
G3D2, Os by u and b in L 2 respectively, integrating by parts over (0,1) x (0, t) with 0 < t < T, 
and using (11.41) i and Lemma l34l we deduce 


J 0 i ^ pu2 + \ h2 + ( x ’t) dx + J o ( a IKI|| 2 + i'||Mia) {x,s)ds 
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Qpo«o + 7^0 + ( x)dx + J v(b 2 (s)b x (l,s) - bi(s)b x {tt, s)) ds 


t r-l 


<C + C I \b(x,t)\dx + C I I (I&2II&I + l^ilH + MH) dxds 
Jo Jo Jo 

<C + C\\b{t)\\ L 2 + C [ (1 + IMlioo) ||6|| L 2ds 
Jo 

< c +-\\b(t )\\ 2 L 2 + c J (1 + II^IIl 2 ) II^IIl 2 ^ 5 

< C +-\\b(t )\\ 2 L 2 + C \\bW\ 2 ds + — J \\u x \\ 2 L2 ds, 

where we have also used Holder’s, Sobolev’s and Cauchy-Schwarz’s inequalities. This, together 
with Gronwall’s inequality, immediately gives (13.61) . □ 

The upper boundedness of the density for the problem (11.41) . (11.101) can be obtained in the 
same manner as in Lemma 12.21 and we have 

Lemma 3.3 Let ( p,u,b ) be a smooth solution of d 1.4j ), hi. 1(A) on [0,1] x [0,T). Then, 

0 <p(x,t)<C, V (x,t) € [0,1] x [0,T). (3.7) 

Next, we aim to estimate first-order derivatives of the velocity and derive the higher inte- 
grability of the magnetic field stated in Lemma 12.51 We remark that it is more complicated to 
achieve this due to the additional boundary effects induced by the magnetic diffusion term vb xx . 
So, to be proceeded, as in Section 2, let “' ” be the material derivative f = ft + uf x and define 


F = ( A u x - P(p) - ^b 2 ) (x, t ) with F x = pit. 


(3.8) 


where ( p,u,b ) is the solution of (11.40 . (11.101) . 

Based on a full use of the “ effective viscous flux” F and subtle boundary analysis, we can 
show the following estimate in a way similar to the proof of Lemma 12.51 

Lemma 3.4 Let ( p,u,b ) be a smooth solution of d 1.4\ ), hi. 1(A) on [0,1] x [0, T). Then, 
sup (|Mt)|||2 + ||6(i)||l 4 + v\\bz(t)\\ 2 L2 ) 


0 <t<T 


+ f (||/3 1 ^ , w|Il 2 + ||&||f,6 + Hl^zlli 2 + V 1 || b xx || £2^ dt < C. 


(3.9) 


Proof. First, similarly to the derivation of (12.111) . (|2.13D and (12.151) . we deduce from () 1.4 1) i and 
(HID 3 that 


A d f 1 2 , [. 2 , 

—r / u x dx + / pu dx 

2 dt ,/n 1 ,/n 


P(p)x(u t + uu x )dx - ^ [ ( b 2 ) x (u t + uu x )dx - ^ [ u x dx 
) 2 Jo 2 j 0 


=sr( pw+ fi“^+ 


r 0 


7 P(p)u 2 x + b 2 u 2 x - ubb xx u x - -u x ) dx 
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(3.10) 


< 


< 


dt 

d_ 

dt 


P{P) + y ) u X ‘ 


v 2, f 1 f 1 

dx + — b 2 x dx + C / (|lire| 3 + ft 6 ) dx + C 

1 Jo Jo 

(^P(p) + — ^ u x dx + — J b 2 xx dx + C J (|i ^| 3 + b G ) dx + C, 


where we have also used (|3.7I) . (13.81) and Cauchy-Schwarz’s inequality. 
Secondly, analogously to the proof of (12.171) . we have 


— [ b A dx + —r f b 6 dx + 12^ [ b 2 b 2 dx 

dt J o 2A J 0 J o 

= j {F + P{p))b 4 dx + [bl(t)b x (l,t) - bl(t)b x (0,t)] 

<C + ^- [ b 6 dx + C [ \F\ 3 dx + Ri{t), 

2A J o Jo 


and consequently, 


d_ 

dt 


f b 4 dx + — f b 6 dx + 12^ f b 2 b 2 x dx < C + C f \F\ 3 dx + Ri(t), (3.11) 

Jo A Jo Jo JO 


where Ri(t) denotes the boundary term: 

Ri(t) = 4z/ [&^(f)6 x (l,f) - bl(t)b x (0,t)] . 

Moreover, multiplying (I1.4ID by 2 vb xx in L 2 and integrating by parts, we use the non-slip 
boundary conditions ul^o.i = 0 to deduce 


b^dx + 2^ 2 / b 2 x dx 


•- r 

dt J 0 

= 2v ( (u x b + ub x ) b xx dx + 2v [b' 2 (t)b x (l,t) - b[ {t)b x (0, t)] 

Jo 

= v / (2 u x bb xx - u x ^) dx + 2z/ [6' 2 (t)6 x (l, t) - 6 / L (t)6 x (0, t)] . 


(3.12) 


Due to non-negativity of P(p ) and 6 2 , it follows from (13.81) that 


- / u x b x dx = —- 


^ + ^(j>) + y )^x<- 


-f 1 

A./o 


Fb x dx , 


which, inserted into (13.121) and combined with (13.71) . (13.81) and Cauchy-Schwarz’s inequality, 
yields 


dt 


b x dx + 2u 2 / b 2 . x dx 


f 

,2 />! 


< 2zz / u T bb TX dx — 


Fb 2 dx + R 2 (t ) 


< 


f 1 b 2 xx dx + C f 1 (ul + b G + u\F\b 2 x ) dx + R 2 (t) 
Jo Jo 
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< 


,2 r 1 /•! 

— / b^ds + C + C ( F 3 + b 6 + is\F\b 2 r ) dx + i? 2 (t), 

2 Jo Jo 


where R 2 (t) is the boundary term: 

#2^) = 2l/ [&' 2 (t)bx(M) - foi(t)6a:(0, *)] . 


Thus, we obtain 


A 

dt 


2 , 3z/ _ 

mdx H—— 


b 2 xx dx <C + C (F 3 + b 6 + i/|F|b£) dx + i? 2 (i). 


(3.13) 


Adding (13.101) and (I3.13P together, we arrive at 

j t J o dx+ f 0 (^ 2 + ^k 2 *) dx 

-JtJo ( P(/ ° ) + t) U * dx + C fj {F 3 + b 6 + v\F\b 2 x )dx + R 2 (t) + C. (3.14) 
Note that it follows from Sobolev’s and Young’s inequalities that for any 0 <t<T, 

\Ri(t)\ + \R 2 (t)\ < Cu\\b x \\ L oo < Cv (||b x || ia + 

< £v 2 \\b X x \\ 2 L 2 + Ce~ x (i'||b x ||| 2 + v 1/2 ^j , (3.15) 


where C > 0 depends on the C 1 -norm of b\(t) and b 2 (t), but is independent of e, v € (0,1). 

Now, multiplying ([3. lift a sufficiently large number K > 0, adding the resulting inequality 
to (13.141) and integrating it over (0,f), by (13.151) with 0 < e < 1 being sufficiently small, we 
deduce in a manner similar to that used in the derivation of (12.20ft that for V v £ (0,1), 


(IM0lll» + IIK0llh + HIW‘)lll0 

+ Jo (^ 1//2 ^^ 2 + l^lli 6 + + ^ll&xzllia) ds 

<C + C f \\Ff L3 ds + Cv [ \\F\\ L -\\b x \\ 2 L 2 ds 

<C + Cv 2 f \\bxt L ,ds + C f {\\F\\ 3 l s + WFW 2 ^) ds 
Jo Jo 

<C + Cv 2 l \\b x \\ 4 L2 ds + i [ \\p 1/2 u\\ 2 L2 ds 
Jo z Jo 

+C f (1 + IWlla) (l + ||«x||i2 + ||b||i 4 ) ds, (3.16) 

Jo 

where (13.61) (13.81) . Sobolev’s, Holder’s and Cauchy-Schwarz’s inequalities have been used to get 


Hollis + imi£~ < c (1 + ||f|| l i) {\\f \\ 2 l2 + \\f\\ l 2 \\f x \\ l2 ) 

< C (1 + KM (l + ||u x ||| 2 + ||&||| 4 ) 

+ C (1 + ||u x ||x,2) (l + ||« x ||l,2 + ||fc|| 2 4) ||p 1 / 2 'd||£ 2 

^ 2 II P 1/2 u\\l 2 + FJ (l + K||| 2 ) (l + K||| 2 + ||b|||4) , 
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which is analogous to (12.221) . As an immediate result of (13.161) . we have 

(IMOIlla + WKVWl* + "IMOIliO 

+ J (\\p 1/2 u \\ 2 L 2 + \\b\\% + v\\bb x \\ 2 L2 + v 2 \\b xx \\ 2 L 2 ^ ds 

< C + Cv 2 f \\b x f Li ds + C f (1 + Hu-eII^) (l + 11 u x ||^2 + H&ll^) ds. 

Jo Jo 

The above inequality, combined with Gronwall’s inequality and the fact that ||m x || 2 2 +^11^x11^2 € 
L 1 (0,T) due to (13.61) . leads to the desired estimate (13.91) . □ 

By virtue of (13.91) . the upper boundedness of the magnetic field and the lower boundedness 
of the density can be obtained in a manner similar to that used in the proof of (12.231) . 

Lemma 3.5 Let ( p,u,b ) be a smooth solution of \l.f\) , on [0,1] x [0, T). Then, 

sup (||6(i)||z,°° + ||/9~ 1 0)||l°°) < C. (3.17) 

0 <t<T 

Proof. Similarly to the proof of (12.241) . we multiply (11.41) -; by 2nb 2n ~ 1 with 1 < n € N and 
integrate by parts over (0,1) to infer, keeping in mind the non-negativity of P(p) and b 2 , that 
there exists a positive constant C , independent of n and u, such that 

— [ b 2n (x, t)dx + 2n(2n — l)u f b 2n ~ 2 b 2 dx 
dt Jo Jo 

= — (2n — 1) f u x b 2n dx + 2nu [62 n_1 (t)6 x (l, t) — 6 2n_1 (t)6 x (0,t)] 

Jo 

< - —— J ^ F + P(p ) + — ^ b 2n dx + CnC 2n v\\b x \\L™ 

Or, — 1 r 1 

< --— / Fb 2n dx + CnC 2n v\\b x \\L°° 

x Jo 

< Cnll-FH^ao||6|| 2 ”„ + CnC 2n u (||6 x||l 2 + IIM^IIM^ 2 ) 

< Cn\\F\\ Laa \\b\\ 2 ^ + CnC 2n (l + v 2 \\b xx \\ 2 L2 ) , 

where we have used (13.91) . Sobolev’s and Cauchy-Schwarz’s inequalities. Noting that urn 
together with (13.81) and (13.91) . gives 

||7 11ioo < \\F\\ 2 L 2 + 11 F 11Li-21| F x 11l 2 < C (1 + II-FxIIl 2 ) 

< c (l + \\p 1/2 u\\ L 2 ^j £L 2 (0,T), (3.18) 

thus we have by (13.91) that 

\\b(t)\\ 2 jfl n < CnC 2n exp i^Cn J \\F\\ L °°dt^ J (l + v 2 \\b xx \\ 2 L2 dt) 

< CnC 2n exp {Cn} , 

from which we conclude that ||&(£)||l°° < C for all t € [0, T) by raising to the power l/(2n) to 
both sides and letting n —» 00. With this at hand, the lower boundedness of the density can be 
shown exactly in the same way as in Lemma 12.61 The proof of (13.171) is complete. □ 
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The following refined estimates of density and magnetic field play an important role in the 
analysis of vanishing resistivity limit, and will be also used to deal with the thickness of magnetic 
boundary layers. 

Lemma 3.6 For any 0 < v < 1. there exists a positive constant C, independent of v, such that 

v l/2 sup (\\b x (t )\\ 2 L 2 + \\px(t)\\ 2 L2 ) + [ (v 3/2 \\b xx \\ 2 L 2 + v 1/2 \\bb x \\ 2 L2 ) dt < C. (3.19) 
0<t<T Jo ' ' 

Proof. Differentiating (11.411 with respect to x, multiplying the resulting equation by 2 p x in L 2 
and integrating by parts, by (13.71) . (13.81) . (|3.17l) and (13.1811 we deduce 


4 [ pl dx = “ 3 [ u x p 2 x dx - 2 f pp x u xx dx 
Jo Jo Jo 

= -3 J u x p 2 dx - Y [ (F x + P x + bb x ) pp x dx 


A _ 

< c (1 + ||||L°°) 11Pa: 11^,2 + C (11^x11^2 + ||Px||l2 + ||&x||l 2 ) IIPzIIl 2 

< C (1 + H^xIIl 2 ) (|| Px II/, 2 + H&xll^) + C (l + ||i7r||| 2 ) , (3.20) 

where we have used the following estimate which follows from (13.71) . (13.81) . (13.171) and (13.181) . 

Il’iprHz. 00 < C (IlFUioo + ||P(p)||z,°° + ||&|||oo) < C (1 + H^xIIl 2 ) ■ (3-21) 

Using the non-slip boundary conditions and integrating by parts, we get from (13.121) that 

d 71 71 71 


dt 


bfdx + 2u z / b xx dx = 2u / (u x b + ub x )b xx dx + R 2 {t) 


= — 2v I u xx bb x dx — 3 z2 / u x b x dx + R 2 (t) + R%(t), 

Jo Jo 

where R 2 (t) is the same as in ()3.131) and i? 3 (i) is also a boundary term denoted by 

i? 3 (f) = 2z2 \b 2 {t){u x b x )(l,t) - bi(t)(u x b x )(0,t)]. 

We are now in a position to bound the terms on the right-hand side of (13.2211 . First, recalling 
the definition of F in (13.811 . using (|3.7j) and (13.171) . we have 

— 2 u / u XT bb T dx = — 


(3.22) 


10 


2v f 

dx = -^- (F x + P(p) x + bb x ) bb x dx 
A Jo 

< Cv(\\f x \\„ + HpxIIl 2 + IIMl 2 ) IIMIzp 

< Cv {\\F x f L2 + 11 Pa: 1)^2 + 11 &X 11 Z.2 ) , 


(3.23) 


and by (13.211) we find 


-3z2 / u x b 2 x dx < Cv\\u x \\ L °°\\b x \\ 2 L 2 < Cv(l+ \\F x \\ L 2 )\\b x \\ 2 L2 . 
Jo 

Choosing 0 < e < 1 small enough, we infer from (13.1511 that 


(3.24) 


R2(t) < ^“ll^xxll^ 2 Cu\\b x \V L 2 + Cu 1/2 . 


(3.25) 
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In order to deal with Rs(t), we observe with the help of (13.211) that 


^x(0,t)| + K(l,£)| < CHuxIIloo < C(1 + \\F x \\ l2 ) , 


and hence, 


R3{t) <Cv( 1 + ||i^||z,2) (||^l ||c([0,T]) + 11^2 ||c([0,T])) HMli 00 

< Cu (1 + H-FxIIl 2 ) (\\b x \\ L 2 + \\b x \\/ 2 11 b xx 11 jji ^ 

2 

— -y ll&xxllla + Cz'HMIl 2 + C (v + u 1 / 2 ^ (l + \\F x \\ 2 l2 ) . (3.26) 

Substituting (13.231) (13.261) into (13.221) . we conclude that for any u € (0,1), 

^IIMi 2 + { l ' 2 \\bxx\\ 2 L2 + HI^MIl 2 ) 

< Cu (||&x|li2 + 11 Px |j^2) (1 + \\F X \\ L .) + Cu 1 / 2 (l + ||iy| 2 ) • (3.27) 

Now, multiplying (13.201) by u and adding the resulting inequality to f|3.27[) . we obtain that 
for any 0 < u < 1, 

(ll^lll 2 + ll^xllia) + u ~\\b X x\\L2 

< Cu (11 Px 11 J;2 + H&xH^) (1 + ||-Fx||l 2 ) + Cv 1 / 2 (l + \\F x \\ 2 l2 ) , 
which, combined with (13.91) . (13.181) and Gronwall’s inequality, implies 

u sup {\\b x \\ 2 L 2 + \\P{p) x \\ 2 L 2 ) + [ (v 2 \\b xx \\ 2 L 2 + u\\bb x \\ 2 L2 )dt 
o <t<T Jo 

< Cu 1 / 2 exp U (1 + II-^x||l 2 ) dt\ J (l + H-Pxll^) dt 

< C'j/ 1 / 2 exp |y (1 + || / o 1 / 2 h|| L 2)dt| J (1 + ||p 1 / 2 h|| 2 2 )dt 

< Cu 1 / 2 , 

which immediately leads to the desired estimate ()3. 19jl . □ 

With the estimates established in Lemmas I2.1H2.81 and 13.2113.61 at hand, we are now ready 
to prove the second part of Theorem 11.21 

Proof of (ii)-Theorem \1.$\ Let (p u , u v , b v ) and ( p,u,b ) be the solutions to the initial¬ 
boundary value problems (11.41) . (11.101) and (11.5(1 . (11.61) . respectively. Then, by a straightforward 
calculation we find that (p v — p,u v — u, b u — b ) satisfies 

(p U - p)t + (p U -PX + P K - n) x + (pT - p) x u" + p x K -u) = 0, (3.28) 

p u (u u - u) t + p v u u (u u - u) x - A ( u v - u) xx = - (p v - p) u t - (p u - p) uu x 

- p" (u- -n)u x - [P(p") - P(p)\ x - i [(6“) 2 - f> 2 ] T , (3.29) 
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and 


(3.30) 


b u 

J XX’ 


(b u - b) t + < ( b v - b) + b (u v - u) x + u v {b» - b) x + (u v -u)b x = ub v x 
First, multiplying (13.3011 by 2(b u — b) in L 2 and integrating by parts, we obtain 

d ^ 


dt 


W~bf L , = - Uv-bf 

Jo 


— b) A dx — 2 / b ( v? — u) x (b v — b) dx 

/0 Jo 

—2 f ( u v — u ) b x ( b v — b) dx + 2v f b xx (b v — b) dx 

Jo Jo 


±h + I 2 + h + h. 


(3.31) 


On one hand, using (12.231) . (12.271) . (|3.18[) . (13.211) and Sobolev’s inequality, we can control 
the first three terms on the right-hand side of (13.311) as follows. 


E^< CIKIU-IK - b \\ 2 L 2 + C\\b\\ LO o\\(u v - u) x \\ L 2 \\b v - b \\ L 2 
i— 1 

+ IIMl 2 II w1/ — u 11 L °°\\b v ~ b\\L 2 

< C (l + II yfff'd \\ L i) || b u — 6111a + 1 ■\\{u u — u) x \\ L 2 \\b u — b \\ L 2 

< e\\{u v — u)x \\ 2 L 2 + Ce 1 (l + Hv^^IIl 2 ) 11^ — 6 |Il 2 - (3.32) 

On the other hand, by (I2.23p . (12.271) and Young’s inequality we deduce 

I A = -2v [ b u x (b u x -b x )dx + 2v{b u x (b u -b){l,t)-b^b"-b)(0,t)] 

Jo 

< —^WKWl 2 + ^II^xIIz^IIMIl 2 + Cv (||6||l°° + ||( 6 i, 62 )||c([o,t])) ll&®IU°° 

< — ^II^xIIl 2 +Cv + Cu ^||6^|| L 2 + ||6£||^ 2 ||6^||J/ 2 ) 

<~^K\\b+ C ^\\^\\b+ Cl/1/2 - (3-33) 

Putting (13.321) and (13.331) into (13.311) . we obtain (0 < e < 1) 

^11^ — &llia + v\\K\\l 2 < C (u 1 ^ 2 + A\Kx\\l 2 ) + e|| (u u - u ) x ||| 2 

+Ce 1 (l + IWW^Wl 2 ) ||6 V — 6|| 2 2 . (3.34) 

Similarly, multiplying (13.281) by 2 (p u — p) in L 2 , integrating by parts, and using Lemmas 
I2.2H2.81 (|3.18l) . (13.211) and Sobolev’s inequality, we deduce 

j t \\p" - p\lb < c|K||i»||/ - p \\ 2 L 2 + cn^iii- \\(u- -uUvW - pIU. 

+ II/L-||l 2 \\ uU ~ u \\l°° || P U ~ P\\l 2 

< c (i + 1| yfifv! Hz, 2 ) II p v - p\\h + C'llK — u )x\\l 2 \\P U — p\\l 2 

< e\\(u u — u) x \\ 2 L 2 + Ce 1 (l + \\\ r p i7 u u \\L 2 ) \\p u ~ p\\ 2 l,2- (3.35) 

Finally, multiplying (13.291) by 2(u u — u) in L 2 and integrating by parts, by virtue of Lemmas 
I2.2H2.81 (13.71) . (13.171) and Sobolev’s inequality again, we see that 

^11 yfWW - «)|lia + 2A||(u I/ - u) x f L2 
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< C(\\u t \\ L 2 + ||u||l<x.||u x || L 2 ) II P v - p\\ L 2\\u v - u\\ L °° + C'lKIUoollv^K - ii)||| 2 

+C (||P(p") - P(p)\\ L 2 + Wf - b 2 || L2 ) ||K - n) x \\ L2 

< C(\\p v - p\\ L 2 + r - b|| i2 ) IIK - «) s || £a + CKIU-Hv^ 7 ^ - u)||| 2 

< e\\(u u — u) x \\ 2 l2 + Ce 1 (Up" — p|| 2 2 + \\b u — &||| 2 ) 

+c (||«x||l 2 + 11'W-xa;11L 2 ) \\VW( uV ~ u )\\ 2 L2 

< e\\(u u — «) x ||^2 + Ce 1 (|| p v — p||| 2 + \\b u — &||| 2 ) + C\\^fp i7 {u 1 ' — u)||^2- (3.36) 

Thus, putting (I3.34[l . (|3.35|) and (13.3611 together, taking £ > 0 suitably small, we arrive at 

^ (II P u ~ p\\h + W ~ 2 + IIVp^K ~ u )\\l' 2 ) + IIK _ u )x\V L 2 + "WKWl 2 

< C (l + IIa/p^u^IIl 2 ) (II P v — p||| 2 + W — &Hl 2 + II\/p"K ~ «)|| l 2 ) 

+C (1M 2 + ^ 2 ||6 xx ||| 2 ) , 

which, combined with (13.91) . ()3. 191) and Gronwall’s inequality, yields 

sup (Up" - p||| 2 + || b u - b\\ 2 L 2 + || yfp(u v - u)||| 2 ) (t) 

0 <t<T 

+ [ (IIK - «)*lli 2 + vW x \\h) dt ^ °u 1/2 . 

■JO 

This inequality, together with the strictly positive lower bound of density due to (13.171) . justifies 
the vanishing resistivity limit stated in the second part of Theorem 11.21 □ 

4 The thickness of magnetic boundary layers 

This section is concerned with the thickness of magnetic boundary layers. As aforementioned in 
Section 1, to simplify the analysis, we assume that the one-dimensional compressible isentropic 
viscous and non-resistive MHD system (jl.5j) is equipped with the specific initial and boundary 
data as follows: 

(p, u, b)\t=o = (p, 0,0) and u| x =o,i = 0, (4.1) 

where p = const, is a positive constant. Then thanks to the uniqueness result, we conclude that 
the initial-boundary value problem (11.51) . (14.11) has only a trivial solution: 

iP,u,b)(x,t) = (p,0,0), (x, t) € [0,1] x [0, T). (4.2) 

This particularly finishes the proof of the first part of Theorem 11.31 

To prove the second part of Theorem 11.31 we consider an initial-boundary value problem 
for (11.41) with initial and boundary data: 

{p,u,b)\ t =o = (p,0,0), u| x=0 ,i = 0, b(0,t) = bi(t) and 6(1, t) = b 2 (t), (4.3) 

where p = const. > 0 and 61 (t), 62 (t) are the same as in (14.11) and (11.101) . respectively. In the 
following, for simplicity but without any confusion, we denote by (p, u, b ) the solution of the 
problem (11.41) . (14.31) . 

Remark 4.1 It is worth mentioning that all the uniform estimates established in Section 3 also 
hold for the solution (p, u, b ) of the problem O, (El- 
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In order to study the thickness of magnetic boundary layers, we need to utilize the trivial 
solution (14.21) of the problem (11.51) . (14.11) . the specific initial data in (14.31) and the convergence 
rates in (|1.14l) to obtain the following improved estimates. 

Lemma 4.1 Let ( p,u,b ) be a smooth solution of 14 -3[ ) on [0,1] x [0,T). Then there exists 

a positive constant C > 0, independent of n, such that 

sup (11/9 -pills + II&III 2 + IMII 2 ) (t) + [ \\u x \\l 2 dt < Cl/ 1 / 2 , (4.4) 

0 <t<T JO 

sup \\u x (t )\\ 2 L 2 + [ Wp^^iiW^dt < Cv 1 / 2 . (4.5) 

0 <t<T Jo 


Proof. As an immediate consequence of (11.141) and (14.21) . one obtains (14.41) . 

To prove (14.51) . similar to the derivation of (13.10ID . we get from (13.71) and (13.171) that 


A d_ 
2 dt 




pu 2 dx 


= f ( P (^) + \ h *) u * dx + J o (t p (p) u I + b2y2 x ~ 

-fjfj (p(p) ~ P(p) T^b 2 ^ n x dx + C (1 + \\u x \\ L °°) 


A o\ 

ubb xx u x — — u% J dx 

11 ^ 11^2 + C^WbxxW^, 


(4.6) 


where we have also used the non-slip boundary conditions u\ x= op = 0. 
Due to the vanishing initial data given in (|4.3D . we have 



(4.7) 


and by (I3.7jl . (13.171) and (14.41) . we get 


P(p) - P{p) + -b 2 ) u x dx 


Thus, by virtue of (14.71) . (14.81) and ()3.191) we 
any 0 < t < T, 


<C(\\p— p\\ L 2 + ||fo||L2) ||u x || L 2 

< “ll^lli 2 + c (II / 5 _ plli 2 + II&III2) 

<j\\u x \\ 2 L 2 + Cu 1/2 . (4.8) 

deduce after integrating (14.61) over (0, t) that for 


u x(f)\\ 2 L,2 + 



ds < Cu 1 ' 2 + C 



U. 


Ux 


I is ds, 


which, together with (13.181) . (13.211) and Gronwall’s inequality, gives (14.51) . □ 

The analysis of the thickness of magnetic boundary layers relies on the following weighted 
(interior) L 2 -estimate of the density and magnetic field. 
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Lemma 4.2 Let ( p,u,b ) be a smooth solution of \l.f\ l, ( f 4-3\) on [0,1] x [0, T). Then there exists 
a positive constant C > 0, independent of v, such that 

sup [ £(x) (\p x \ 2 + \b x \ 2 ) (x,t)dx < Cv l/2 , (4.9) 

0<t<T Jo 


where £(x) = x 2 (l — x) 2 . 

Proof. Differentiating (11.41) ^ with respect to x gives 


b x t — vbxxx = —ub xx — 2u x b x — u xx b, 


which, multiplied by b x £,(x) in L 2 and integrated by parts, yields 

bU(x)dx + vJ^ b 2 xx i{x)dx 

= ^ J b 2 f"(x)dx + ub 2 r f'(x)dx 

3 rl r 1 ^ ^ 

“o / u x t> 2 x £(x)dx- / u xx bb x f(x)dx = ^ I t , 
1 Jo Jo 


(4.10) 


where we have used the fact that £(x) = f'(x) = 0 on the boundaries x = 0,1. 

We are now in a position of estimating each term on the right-hand side of (14.101) . First, it 
readily follows from (13.191) that 


h < Cv\\b x \\ 2 L 2 < Cv 1/2 . 


(4.11) 


Noting that due to the non-slip boundary conditions u\ x= 03 = 0 and (|4.5I) (also cf. (13.91) ), 
one has 

px 

u(x,t) = / u x (x,t)dx < Cx\\u x \\l°° , ViG [0,1] 

Jo 

u(x,t) = — u x (x,t)dx < C(1 — x)\\u x \\l°°, V i G [0,1], 

J X 

so that, by (|3.18D and (|3.21l) . we have 

[ ub 2 [x(l — x) 2 — x 2 (l — x)] dx 

Jo 


h = 


< C\\u x \\l°° / b x f{x)dx 

Jo 

< C (l + ||p 1/2 ii|| L2 ) J b 2 x Z(x)cL 


(4.12) 


Similarly, 


h < C\\u x \\ L o° J h 2 x i{x)dx <C + \\p l/2 u\\ L 2 ^j J b 2 £(x)dx. (4.13) 
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Finally, recalling that 


Uxx — A (Fx T P(p)x “I - bbx') , F x — pii , 
and using (13.71) . (13.17H and Cauchy-Schwarz’s inequality, we infer from (13.81) that 

-A' 1 / (F x + P(p) x + bb x ) bb x £(x)dx 
Jo 


h = -; 


< 


< 


< 


c r (|F X | + \ Px \ + \b x \) \b x Ux)dx 
Jo 

C\\F x \\\ 2 + C f 1 (pi + b 2 x ) t(x)dx 
Jo 

C\\p 1/2 u \\ 2 L 2 +C f ( p 2 x + b 2 x ) Z{x)dx. 

Jo 


(4.14) 


Thus, inserting (14. 11 D — (14.14D into (14.101) . we arrive at 
d n n 


dt So h l^ x ) dx + v J 0 b lxt( x )dx 

< C (v 1/2 + \\p 1/2 u\\ 2 L2 ^j +C (l + \\p 1/2 u\\ L 2 ^j j ( pi + bl) £(x)dx . (4.15) 

In a similar manner, differentiating (11.41) i with respect to x, multiplying it by p x £(x) in L 2 , 
and integrating by parts, we deduce 


1 d_ 

2 dt 


[ plC{x)d: 

Jo 


X = 


< 


1 3 r l r 1 

up 2 x i\x)dx - - u x p x £(x)dx - / u xx pp x £(x)da 
z Jo Jo 

C||^||| 2 +C(1 + |K|| L ~) [ (, p 2 x + bl)Z(x)da 

Jo 


< ( , ||p 1 ^ 2 h ||^2 + C ^1 + ||p 1 ^ 2 h|| I/ 2 ^) J (p 2 + b 2 ) £(x)dx, 

which, combined with ()4. 151) . shows that 

j t J {pI + b 2 x ) £(x)dx + v J b 2 xx £(x)dx 
< c (y l/2 + \\p 1/2 d\\ 2 L 2 ^ +c (l + ||pi/ 2 ft|| L 2 ) J (pl + bl) t(x)dx. (4.16) 


Therefore, by virtue of (14.51) . we deduce from (14.161) that 
/ £(*) (pi + h l ) (x, t)dx + u [ 

0 <t<T 



0 JO 


sup [ f (x) (pi + 6 2 ) {x, t)dx + 

<t<T Jo 

<C^ 2 + J \\p 1/2 u\\ 2 L 2dtj exp^C J (^1 + ||p 1/2 tt|| L 2^ dtj 


£( x)b xx (x,t)dxdt 
T 


< Cv 1 / 2 . 

The proof of (14.91) is complete. 


□ 
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With the help of (14.41) and (14.91) . we can estimate the thickness of magnetic boundary layers. 
Proof of fii)-Theorem \ 1.31 On one hand, by (14.91) we observe that for any 5 G (0,1/2), 


<5 2 [ (Pi + b 2 x ) (x, t)dx 

4 6 

rl/2 pi-8 

= b 2 (Pl + b%) (x, t)dx + 5 2 ( p 2 x + h 2 x ) (x, t)dx 

J8 J 1/2 

I* 1/2 i*l—8 

< X 2 (pi + h 2 x ) (x, t)dx + / (1 - x) 2 (pi + hi) (x, t)<k 

J 5 J 1/2 

- 1/2 


< 4 


< 4 


/•I/Z l-L—O 

■ / a;2 (l - x ) 2 (p 2 + €) (*, t)dx + 4 / x 2 (l - x) 2 (p 2 + hi) (x, t)da 

45 41/2 

J x 2 (l — x) 2 (pi + h 2 x ) (x, t)dx < Cv 1 / 2 , 

from which it follows that for any 5 G (0,1/2), 

\\Px(t)\\L 2 (n 5 ) + \\bx(t)\\L 2 (n s ) — C5 V 1 / 4 , VtG[0,T), 


(4.17) 


where — (<5,1 — <5). Thus, in view of (14.41) and (14.171) . we have by Sobolev’s inequality that 
for any t G [0, T), 

\\(p-m)\\ 2 c (ns) + \m\\ 2 c( n s) 

- C (ll(p-p)( i )l| 2 2(n 5 ) + IIKOH^n*)) 

+C (IK/ 3 _ P)( i )llL 2 (o i )ll7 , .T(i)l|L 2 (o i ) + H^( i )llL 2 (n 5 )ll^(^)llL 2 (n 5 )) 

< C (ll(p — p){t) lli 2 (n) + IK(/)lll 2 (Q)) 

+C* (IK/ 3 _ P) (/) Hl 2 (o) ll/ 3 *(^) lli/ 2 (Oa) + IK(^) IIl 2 (h) ||^x(^) 11X/ 2 )) 

< CV 1/2 + Cv 1 ' 2 ^ 1 < Cv 1/2 5~ l -> 0 as v -> 0, (4.18) 


provided 6 = 5{v) satisfies 


5(r^) —>- 0 and 


6(y)_ 

1/1/2 


oo, 


On the other hand, it is easy to see that 


as v —>• 0 . 


(ll(P P)(^)llicx,(o,T;C(n)) ll^(^llioo( 0T;C /(Q))) >0) 

provided the boundary data bi(t),b- 2 (t) are not identically zero (i.e., b\(t), 62 (f) / 0). This, 
together with (|4.18l) , proves the second part of Theorem 11.31 □ 
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